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Abstract. On the way to solid-state quantum computing, overcoming decoherence 
is the central issue. In this contribution, we discuss the modeling of decoherence of 
a superonducting flux qubit coupled to dissipative electronic circuitry. We discuss 
its impact on single qubit decoherence rates and on the performance of two-qubit 
gates. These results can be used for designing decoherence-optimal setups. 

Quantum computation is one of the central interdisciplinary research themes 
in present-day physics [1]. It promises a detailed understanding of the often 
counterintuitive predictions of basic quantum mechanics as well as a quali- 
tative speedup of certain hard computational problems. A generic, although 
not necessarily exclusive, set of criteria for building quantum computers has 
been put forward by DiVincenzo [2]. The experimental realization of quan- 
tum bits has been pioneered in atomic physics, optics and NMR. There, the 
approach is taken to use microscopic degrees of freedom which are well iso- 
lated and can be kept quantum coherent over long times. Efficient controls 
are attached to these degrees of freedom. Even though these approaches are 
immensely succesful demonstrating elementary operations, it is not evident 
how they can be scaled up to macroscopic computers. 

Solid-state systems on the other hand have proven to be scalable in 
present-day classical computers. Several proposals for solid-state based quan- 
tum computers have been put forward, many of them in the context of su- 
perconductors [3]. As solid-state systems contain a macroscopic number of 
degrees of freedom, they are very sensitive to decoherence. Mastering and 
optimizing this decoherence is a formidable task and requires deep under- 
standing of the physical system under investigation. Recent experimental 
success [4,5] suggests that this task can in principle be performed. 

In this contribution, we are going to study decoherence of superconducting 
qubits coupled to an electromagnetic environment which produces Johnson- 
Nyquist noise. We show, how the decoherence properties can be engineered 
by carefully designing the environmental impedance. We will discuss how the 
decoherence affects the performance of a CNOT operation. 



2 



F.K. Wilhclm et al. 



1 Superconducting Flux Qubits 

Superconducting qubits [3-6] are very well suited for the task of solid-state 
quantum computation, because two of the most obvious decoherence sources 
in solid-state systems are supressed: Quasiparticle excitations experience an 
energy gap and phonons are frozen out at low temperatures [7]. The com- 
putational Hilbert space is engineered using Josephson tunnel junctions that 
are characterized by two competing energy scales: The Josephson coupling 
of a junction with critical current I c , Ej = I c ^q/2tt, and the charging en- 
ergy E c h = 2e 2 /Cj of a single Cooper pair on the geometric capacitance 
Cj of the junction. Here <?o = h/2e is the superconducting flux quantum. 
There is a variety of qubit proposals classified by the ratio of this ener- 
gies. Whereas another contribution in this volume [8] focuses on the case 
of charge qubits, E c h > Ej, this contribution is motivated by flux qubit 
physics, Ej > E^. However, most of the discussion has its counterpart in 
other superconducting setups as well. Specifically, we discuss a three junction 
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Fig. 1. Experimental setup for measurements on a flux qubit. The qubit (center) 
is a superconducting loop that contains three Josephson junctions. It is induc- 
tively coupled to a DC-SQUID (a), and superconducting control lines for applying 
magnetic fields at microwave frequencies (b) and static magnetic fields (c). The 
DC-SQUID is realized with an on-chip shunt circuit with impedance Z(ui). The 
circuits a)-c) are connected to filtering and electronics (not drawn) 



qubit [6,9], a micrometer-sized low- inductance superconducting loop contain- 
ing three Josephson tunnel junctions (Fig. 1). By applying an external flux 
<P q a persistent supercurrent can be induced in the loop. For values where <P q 
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is close to a half-integer number of flux quanta, two states with persistent 
currents of opposite sign are nearly degenerate but separated by an energy 
barrier. We will assume here that the system is operated near <I> q — ^<Po. The 
persistent currents in the classically stable states have here a magnitude I p . 
Tunneling through the barrier causes a coupling between the two states, and 
at low energies the loop can be described by a Hamiltonian of a two state 
system [6,9], 

e A 

H q = -a z + — a x , (1) 

where a z and a x are Pauli matrices. The two eigenvectors of a z correspond 
to states that have a left or a right circulating current and will be denoted as 
\L) and \R). The energy bias e = 2I p (<P q — ^<P ) is controlled by the externally 
applied field <P q . We follow [10] and define A as the tunnel splitting at $ q = 
i^Oj such that A — 2W with W the tunnel coupling between the persistent- 
current states. This system has two energy eigen values A 2 + e 2 , such 
that the level separation v gives v = V A 2 + e 2 . In general A is a function 
of e. However, it varies on the scale of the single junction plasma frequency, 
which is much above the typical energy range at which the qubit is operated, 
such that we can assume A to be constant for the purpose of this paper. 

In the experiments <^ q can be controlled by applying a magnetic field 
with a superconducting coil at a distance from the qubit and for local control 
one can apply currents to superconducting lines, fabricated on-chip in the 
vicinity of the qubit. The qubit 's quantum dynamics will be controlled with 
resonant microwave pulses (i. e. by Rabi oscillations). In recent experiments 
the qubits were operated at e ss 5Z\ ore~0 [4,9]. The numerical values given 
in this paper will concentrate on the former case. At this point, there is a 
good trade-off between a system with significant tunneling, and a system with 
cr z -like eigenstates that can be used for qubit-qubit couplings and measuring 
qubit states [6] . The qubit has a magnetic dipole moment as a result of the 
clockwise or counter-clockwise persistent current The corresponding flux in 
the loop is much smaller than the applied flux <P q , but large enough to be 
detected with a SQUID. This will be used for measuring the qubit states. For 
our two- level system Eq. (1), this means that both manipulation and readout 
couple to <7 Z . Consequently, the Nyquist noise produced by the necessary 
external circuitry will couple in as flux noise and hence couple to a z , giving 
e a small, stochastically time-dependent part Se(t). 

Operation at e w has the advantage that the flux noise leads to less 
variation of v. In the first experiments [4] this has turned out to be crucial 
for observing time-resolved quantum dynamics. Here, the qubit states can be 
measured by incorporating the qubit inside the DC-SQUID loop. While not 
working that out in detail, the methods that we present can also be applied 
for the analysis of this approach. This also applies to the analysis of the 
impact of electric dipole moments, represented by u x . With E c h <C Ej, these 
couple much less to the circuitry and will hence not be discussed here. 
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As the internal baths are well suppressed, the coupling to the electromag- 
netic environment (circuitry, radiation noise) becomes a dominant source of 
decoherence. This is a subtle issue: It is not possible to couple the circuitry 
arbitrarily weakly or seal the experimental setup, because it has to remain 
possible to control the system. One rather has to engineer the electromagnetic 
environment to combine good control with low unwanted back-action. 

Any linear electromagnetic environment can be described by an effective 
impedance Z e g. If the circuit contains Josephson junctions below their crit- 
ical current, they can be included through their kinetic inductance Lkin = 
<1>o/(2itI c cos(f>), where (f> is the average phase drop across the junction. The 
circuitry disturbs the qubit through its Johnson-Nyquist noise, which has 
Gaussian statistics and can thus be described by an effective Spin-Boson 
model [11]. In this model, the properties of the oscillator bath which forms 
the environment are characterized through a spectral function J (to), which 
can be derived from the external impedance. Note, that other nonlinear ele- 
ments such as tunnel junctions which can produce non-Gaussian shot noise 
are generically not covered by oscillator bath models. 

As explained above, the flux noise from an external circuit leads to e = 
eo + Se(t) in Eq. (1). We parametrize the noise Se(t) by its power spectrum 



Thus, from the noise properties calculated by other means one can find J(uj) 
as was explained in Detail in [12]. In this contribution, we would like to out- 
line an alternative approach pioneered by Leggett [13], where J(w) is derived 
from the classical friction induced by the environment. In reality, the com- 
bined system of SQUID and qubit will experience fluctuations arising from 
additional circuit elements at different temperatures, which can be treated 
in a rather straightforward manner. 

2 Decoherence from the Electromagnetic Environment 

2.1 Characterizing the Environment from Classical Friction 

We study a DC-SQUID in an electrical circuit as shown in Fig. 1. It contains 
two Josephson junctions with phase drops denoted by 7i/ 2 - We start by 
looking at the average phase 7 ox = (71 + 72)/2 across the read-out SQUID. 
Analyzing the circuit with Kirchhoff rules, we find the equation of motion 



2Cj-^7e X = -2/ c ,o cos( 7j ) sin 7cx + / bias - / dt'% x (t')Y{t - t'). (3) 



Here, 7m = (71 — 72) /2 is the dynamical variable describing the circulating 
current in the loop which is controlled by the flux, /bias is the bias current 
imposed by the source, Y(ui) = Z~ x {lo) is the admittance in parallel to the 
whole SQUID and Y(j) its Fourier transform. The SQUID is described by 



({<Se(t),<fe(0)}) u = h 2 J(w)coth(hu/2k B T). 



(2) 
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the junction critical currents I c ,o which arc assumed to be equal, and their 
capacitances Cj. We now proceed by finding a static solution which sets the 
operation point 7i n / e x,o an d small fluctuations around them, c>7i n / ox . The 
static solution reads /bias = I c ,ctf sin7 ox ,o where 7 Cie ff = 2/ c , cos7i„,o is the 
cfFcctive critical current of the SQUID. Linearizing Eq. 3 around this solution 
and Fourier-transforming, we find that 

27r/ b tan7 in Z cS (uj) 
<57cxM = — ^ S^u) (4) 

where Z e s(u>) = (Z(cj) -1 + 2iu>Cj + (iwLkin) 1 ) 1 is the effective impedance 
of the parallel circuit consisting of the Z(u>), the kinetic inductance of the 
SQUID and the capacitance of its junctions. Neglecting self-inductance of the 
SQUID and the (high-frequency) internal plasma mode, we can straightfor- 
wardly substitute 7; n = 7r<£/<£ and split it into 7^0 = tt^x,s/^o set by the 
externally applied flux <£ x g through the SQUID loop and Sji = ttMsqIq/<P 
where MgQ is the mutual inductance between qubit and the SQUID and 
Iq(<P) is the circulating current in the qubit as a function of the junction 
phases, which assumes values ±J p in the classically stable states. 

In order to analyze the backaction of the SQUID onto the qubit in the 
two-state approximation, Eq. (1), we have to get back to its full, continuous 
description, starting from the classical dynamcis. These are equivalent to a 
particle, whose coordinates are the two independent junction phases in the 
three-junction loop, in a two-dimensional potential 

C($ /2n) 2 (p = -VU(tp, $ x , g + I S M SQ ). (5) 

The details of this equation are explained in [6] . C is the capacitance matrix 
describing the charging of the Josephson junctions in the loop, U (if) con- 
tains the Josephson energies of the junctions as a function of the junction 
phases and Is- is the ciculating current in the SQUID loop. The applied flux 
through the qubit ^ q is split into the flux from the external coil <P x ,q and the 
contribution form the SQUID. Using the above relations we find 

I S M SQ = 6$ cl - 2tt 2 M s 2 q /! tan 2 7in , (6) 

where 8<P C \ ~ Msq/ Cj0 cos7 cx .o sin7; nj0 is the non-fluctuating back-action 
from the SQUID. 

From the two-dimensional problem, we can now restrict ourselves to the 
onc-dimcnsional subspace defined by the preferred tunneling direction [6], 
which is described by an effective phase ip. The potential restricted on this 
direction, Um{tp) has the form of a double well [11,14] with stable minima 
situated at ±(fio- In this way, we can expand Z7iD(y,^ q ) — U(<p 7 <P q ,x) + 
Iq(<p)IqMsq. Approximating the phase-dependence of the circulating current 
as Iq((p) ~ Ipip/ipo where I p the circulating current in one of the stable 
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minima of ip, we end up with the classical equation of motion of the qubit 
including the backaction and the friction induced from the SQUID 



= -d lp Uir>(<P,$x,q 



27r 2 M s 2 Q 4 2 i: 
- 5* c i). 



,tan 7in,0T 



(7) 



From this form, encoded as D{u>)ip(uj) = —dU/dtp we can use the prescription 
given in [13] and identify the spectral function for the continuous, classical 
model as J con t = Iid-D(lu). From there, we can do the two-state approximation 
for the particle in a double well [14] and find J(w) in analogy to [12] 



J(u 



(2tt) 2 /M SQ / F 



^bias * an2 



—J Re{Z cS (cj)}. 



(8) 



2.2 Qubit Dynamics under the Influence of Decoherence 

From J(u), we can analyze the dynamics of the system by studying the 
reduced density matrix, i.e. the density matrix of the full system where the 
details of the environment have been integrated out, by a number of different 
methods. The low damping limit, J{uo)/u <C 1 for all frequencies, is most 
desirable for quantum computation. Thus, the energy-eigenstates of the qubit 
Hamiltonian, Eq. (1), are the appropriate starting point of our discussion. In 
this case, the relaxation rate r r (and relaxation time r r ) are determined 
by the environmental spectral function J(oj) at the frequency of the level 
separation v of the qubit 




where T is the temperature of the bath. The dephasing rate (and dephas- 
ing time r^) is 

r,- V - -£ + ara (£)' ,io, 

with a = lim^^o J(u))/(2ircu). These expressions have been derived in the 
context of NMR [15] and recently been confirmed by a full path-integral 
analysis [10]. In this paper, all rates are calculated for this regime. 

For performing efficient measurement, one can afford to go to the strong 
damping regime. A well-known approach to this problem, the noninteracting 
blip approximation (NIBA) has been derived in Ref. [13]. This approximation 
gives good predictions at degeneracy, e = 0. At low |e| > it contains an 
artifact predicting incoherent dynamics even at weak damping. At high bias, 
e^> A and at strong damping, it becomes asymptotically correct again. We 
will not detail this approach here more, as it has been extensively covered in 
[11,14]. 
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If J(u>) is not smooth but contains strong peaks the situation becomes 
more involved: At some frequencies, J(u>) may fall in the weak and at others 
in the strong damping limit. In some cases, whcrn J(u>) <C w holds at least 
for uj < ft with some ft ^> v/K, this can be treated approximately: one can 
first renormalize A e g through the high-frequency contributions [11] and then 
perform a weak-damping approximation from the fixed-point Hamiltonian. 
This is detailed in Ref. [16]. In the general case, more involved methods such 
as flow equation renormalization [17] have to be used. 

3 Engineering the Measurement Apparatus 

From Eq. 8 we see that engineering the decoherence induced by the mea- 
surement apparatus essentially means engineering Z c g. This includes also 
the contributions due to the measurement apparatus. In this section, we are 
going to outline and compare several options suggested in literature. We as- 
sume a perfect current source that ramps the bias current /bias through the 
SQUID. The fact that the current source is non-ideal, and that the wiring to 
the SQUID chip has an impedance is all modeled by the impedance Z(ui). 
The wiring can be engineered such that for a very wide frequency range 
the impedance Z(w) is on the order of the vacuum impedance, and can be 
modeled by its real part Ri. It typically has a value of 100 ft. 

3.1 R-Shunt 

It has been suggested [18] to overdamp the SQUID by making the shunt cir- 
cuit a simple resistor Z{lo) = R$ with i?s <C \/£kin/2Cj. This is inspired by 
an analogous setup for charge qubits, [3]. Following the parameters given in 
[12], a SQUID with J C;0 = 200nA at <P/<P ^ 0.75 biased at / b ias = 120nA, we 
find L kin ~ 2 • 10" 9 iJ. Together with Cj ~ lfF, this means that the SQUID 
is overdamped if R <C i? m ax = 1.4kfi. Using Eq. 8, we find that this provides 
an Ohmic environment with Drude-cutoff, J(uj) = au)/(l + lo 2 /lo 2 lr ) where 
lu lr = R/L kia and a = {2ir) 2 /h {M SQ I q /$ f I b 2 ias ta,n 2 (n$ /$ Q )LlJR s . Us- 
ing the parameters from Ref. [12], M SQ I q /$o = 0.002, we find aR = 0.08fJ 
and u)lr/R = 8.3GHz/f2. Thus, for our range of parameters (which essen- 
tially correspond to weak coupling between SQUID and qubit), one still has 
low damping of the qubit from the (internally overdamped) environment at 
reasonable shunt resistances down to tens of Ohms. For such a setup, one 
can apply the continuous weak measurement theory as it is outlined e.g. in 
[18]. This way, one can readily describe the readout through measurement of 
Z e g which leaves the system on the superconducting branch. If one desires to 
read out the state by monitoring the voltage at bias currents above the I c , e s, 
our analysis only describes the pre-measurement phase and at least shows 
that the system is hardly disturbed when the current is ramped. 
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Fig. 2. Circuit models for the C-shunted DC-SQUID (a) and the 7?C-shunted 
DC-SQUID (b). The SQUID is modeled as an inductance Lj. A shunt circuit, the 
superconducting capacitor C s h or the -R s h-C s h series, is fabricated on chip very 
close to the SQUID. The noise that couples to the qubit results from Johnson- 
Nyquist voltage noise 8V from the circuit's total impedance Z e a. Z e a is formed 
by a parallel combination of the impedances of the leads Zi, the shunt and the 
SQUID, such that Z~g = 1/Z t + l/(R sh + l/icuC sh ) + 1/icuL.j, with R sh = for 
the circuit (a) 



3.2 Capacitive Shunt 



Next, we consider a large superconducting capacitive shunt (Fig. 2a, as im- 
plemented in Refs. [4,9]). The C shunt only makes the effective mass of the 
SQUID's external phase 7 ex very heavy. The total impedance Z c s(lj) and 
J{lo) are modeled as before, see Fig. 3. As limiting values, we find 



Re{Z eff (u;)} 



Ri ' 

Ri-, 
i 



for u <C uolc 
for ijj ~ ujlc 
for uj » uolc 



(11) 



We can observe that this circuit is a weakly damped LC-oscillator and 
it is clear from (9) and (8) that one should keep its resonance frequency 
w lc = l/V-^jCsh, where Re{Z e s(ui)} has a maximum, away from the qubit's 
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co/271 (Hz) oo/2jc (Hz) 

Fig. 3. A typical Re{Z t (ui)} for the C-shunted SQUID (a) and the .RC-shunted 
SQUID (b), and corresponding J{u>) in (c) and (d) respectively. For comparison, 
the dashed line in (c) shows a simple Ohmic spectrum, J{u>) = aui with exponential 
cut off uJc/2-k = 0.5 GHz and a = 0.00062. The parameters used here are I p — 500 
nA and T — 30 mK. The SQUID with 2I co = 200 nA is operated at / = 0.75 tt and 
current biased at 120 nA, a typical value for switching of the C-shunted circuit 
(the 7?C-shunted circuit switches at higher current values). The mutual inductance 
M = 8 pH (i. e. MI p /$o = 0.002). The shunt is C sh = 30 pF and for the RC shunt 
R sh = 10 fi. The leads are modeled by Ri = 100 fi 



resonance w ros = v/K. This is usually done by chosing wlc "C w rC s- For a 
C-shunted circuit with ulc "C w res , this yields for J(uj w wlc) 

T . , (2tt) 2 /M7 p \ 2 t2 o /tt^\ 1 

The factor l/o> 3 indicates a natural cut-off for J(w), which prevents the ultra- 
violet divergence [11,10] and which in much of the theoretical literature is in- 
troduced by hand. Using Eq. 9, we can directly analyze mixing times r r vs oj ics 
for typical sample parameters (here calculated with the non-approximated 
version of Rje{Z t (u))}) , see Ref. [12] for details. The mixing rate is then r r w 
(2ttA/H) 2 w -» (MI p /$ f tan 2 (7r^/<l> )(2?iC s \i? ; )- 1 coth (hcj res /2k B T) . 
With the C-shunted circuit it seems possible to get t t values that are very 
long. They are compatible with the ramp times of the SQUID, but too slow 
for fast repetition rates. For the parameters used here they are in the range of 
15 /is. While this value is close to the desired order of magnitude, one has to 
be aware of the fact that at these high switching current values the lineariza- 
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tion of the junction as a kinetic inductor may underestimate the actual noise. 
In that regime, phase diffusion between different minima of the washboard 
potential also becomes relevant and changes the noise properties [19,20]. 

3.3 RC-Shunt 

As an alternative we will consider a shunt that is a series combination of a 
capacitor and a resistor (Fig. 2b) (i?C-shunted SQUID). The RC shunt also 
adds damping at the plasma frequency of the SQUID, which is needed for 
realizing a high resolution of the SQUID readout (i. e. for narrow switching- 
current histograms) [19]. The total impedance Z t {to) of the two measurement 
circuits are modeled as in Fig. 2. For the circuit with the RC shunt 



Re{Z t (w)} 



The difference mainly concerns frequencies ui > lolc, where the C-shunted 
circuit has a stronger cutoff in Re{Z ff (w)}, and thereby a relaxation rate, 
that is several orders lower than for the i?C-shunted circuit. Given the values 
of J(u>) from Fig. 3 one can directly see from the values of that an RC- 
shunted circuit with otherwise similar parameters yields at u res /2-K = 10 
GHz relaxation times that are about four orders of magnitude shorter. 

4 Coupled Qubits 

So far, we have applied our modeling only to single qubits. In order to study 
entanglement in a controlled way and to eventually perform quantum algo- 
rithms, this has to be extended to coupled qubits. 

4.1 Hamiltonian 

There is a number of ways how to couple two solid-state qubits in a way which 
permits universal quantum compuation. If the qubit states are given through 
real spins, one typically obtains a Heisenberg-type exchange coupling. For 
other qubits, the three components of the pseudo-spin typically correspond 
to physically completely distinct variables. In our case, a z corresponds to the 
flux through the loop whereas <r x / y are charges. Consequently, one usually 

finds Ising-type couplings. The case of a [ y ] ® af ] coupling, i.e. coupling by 
a component which is orthogonal to all possible single-qubit Hamiltonians, 
has been extensively studied [21,22], because this type is straightforwardly 
realized as a tunable coupling of charge qubits [3] . We study the generic case 
of coupling the "natural" variables of the pseudospin to each other, which 



Ri 



for u) <C uo LC 



<R U for ui = u LC « fl-W (13) 
Rl/ /Rah, for U) = UJLC > 



RshC 

Ril/Rsh, forw>w LC 
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can be realized in flux qubits using a switchable superconducting transformer 
[6,23], but has also been experimentally utilized for coupling charge qubits 
by fixed capacitive interaction [24] . 

We model the Hamiltonian of a system of two qubits, coupled via Ising- 
type coupling. Each of the two qubits is described by the Hamiltonian Eq. 
(1). The coupling between the qubits is described by H qq = —(K/2)a^ <g> 

(2) 

ai that represents e.g. inductive interaction. Thus, the complete two-qubit 
Hamiltonian in the absence of a dissipative environment reads 

Kb = -\ £ + At*?) \Ka?>vV>. (14) 

i=l,2 

For two qubits, there are several ways to couple to the environment: Both 
qubits may couple to a common bath such as picked up by coupling elements 
[6]. Local readout and control electronics coupling to individual qubits [6] 
can be described as coupling to two uncorrelated baths. In analogy to the 
procedure described above, one can determine the spectral functions of these 
baths by investigating the corresponding impedances. 

In the case of two uncorrelated baths, the full Hamiltonian reads 

K = K + ^ * W + K + H B2 , (15) 

i=l,2 

_X"W = CYl<v ^v x v are collective coordinates of the bath. In the case of two 
qubits coupling to one common bath we model our two qubit system in a 
similar way with the Hamiltonian 

Kb = Kb + \ (*P + ^ ] ) X + H B 
where X is a collective bath coordinate similar to above. 

4.2 Rates 

We can derive formulae for relaxation and dephasing rates similar to Eqs. (9) 
and (10). Our Hilbert space is now four-dimensional. We label the eigenstates 
as \E1) . . . \E4). We chose \E1) to be the singlet state (|T|) - ||t)) /\/2, which 
is always an eigenstate [25] whereas \E2) . . . \E4) are the energy eigenstates 
in the triplet subspace, which are typically not the eigenstates of 
As we have 4 levels, we have 6 independent possible quantum coherent os- 
cillations, each of which has its own dephasing rate, as well as 4 relaxation 
channels, one of which has a vanishing rate indicating the existence of a stable 
thermal equilibrium point. The expressions for the rates, although of simi- 
lar form as in Eqs. (9) and (10) are rather involved and are shown in [25]. 
Figure 4 displays the dependence of typical dephasing rates and the sum of 
all relaxation rates on temperature for the case A = e = K = hi>s with 



(16) 
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i>s = 1GHz. The rates are of the same magnitude for the case of one common 
bath and two distinct baths. If the temperature is increased above the roll 
off point set by the intrinsic energy scales, T s = (h/ks)^ = 4.8 • 10~ 2 K, 
where E s = 1GHz, the increase of the dephasing and relaxation rates follows 
a linear dependence, indicating that the environmental fluctuations are pre- 
dominantly thermal. As a notable exception, in the case of one common bath 




Fig. 4. Log- log plot of the temperature dependence of the sum of the four relaxation 
rates and selected dephasing rates. Qubit parameters K, e and r\ are all set to E s and 
the bath is assumed to be Ohmic a — 10 -3 . The upper panel shows the case of one 
common bath, the lower panel the case of two distinct baths. At the characteristic 
temperature of approximately 0.1 • T s the rates increase very steeply 



the dephasing rates r v21 = r vi2 go to zero when the temperature is decreased 
while all other rates saturate for T — > 0. This can be understood as follows: 
the singlet state \E1) is left invariant by the Hamiltonian of coupled qubits 
in a common bath, Eq. (16), i.e. it is an energy eigenstate left unaffected 
by the environment. Superpositions of the singlet with another eigenstate 
are usually still unstable, because the other eigenstate generally suffers from 
decoherence. However, the lowest-energy state of the triplet subspace \E2) 
cannot decay by spontaneous emission and flip-less dephasing vanishes at 
T = 0, hence the dephasing rate between eigenstates \E1) and \E2) vanishes 
at low temperatures, see Fig. 4. As shown in [25], there can be more "pro- 
tected" transitions of this kind if the qubit parameters are adjusted such that 
the symmetry between the unperturbed qubit and the coupling to the bath 
is even higher, e.g. at the working point for a CPHASE operation. 
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4.3 Gate Performance 

The rates derived in the previous section are numerous and do strongly de- 
pendend on the tunable parameters of the qubit. Thus, they do not yet allow 
a full assesment of the performance as a quantum logic element. A quan- 
titative measure of how well a two-qubit setup performs a quantum logic 
gate operation are the gate quality factors introduced in [26]: the fidelity, 
purity, quantum degree and entanglement capability. These factors charac- 
terize the density matrices obtained after attempting to perform the gate 
operation in a hostile environment, starting from all possible initial condi- 
tions p(0) = \&l n ) (&i n \ - To form all possible initial density matrices needed to 
calculate the gate quality factors, we use the 16 unentangled product states 
\&i n ), j = 1,. .. ,16 defined [22] according to \V a ) 1 \%) 2 , {a,b= 1, . . . , 4), with 
= |0), |# 2 ) = |1>, |# 3 ) = (1A/2)(|0> + |1», and |# 4 ) = (l/>/2)(|0> +i |1)). 
They form one possible basis set for the superoperator v G which describes the 
open system dynamics such that p(to) — ^gp(O) [22,26]. The CNOT gate is 
implemented using rectangular DC pulses and describing dissipation through 
the Bloch-Redfield equation as described in [3,25]. 

The fidelity is defined as T = (l/16)£)]=i (&L\ U G Pg U g\&L) • The fi- 
delity is a measure of how well a quantum logic operation was performed. 
Clearly, the fidelity for the ideal quantum gate operation is equal to 1. The 
second quantifier is the purity V = (1/16) Y^Li tr (Pg) 2 > which should be 
1 in a pure and 1/4 in a fully mixed state. The purity characterizes the effects 
of decoherence. The quantum degree measures nonlocality. It is defined as the 
maximum overlap of the resulting density matrix after the quantum gate op- 
eration with the maximally entangled Bell-states Q = maxj.fe {^meiPG^mei '■ 
For an ideal entangling operation, e.g. the CNOT gate, the quantum degree 
should be 1. It has been shown [27] that all density operators that have an 
overlap with a maximally entangled state that is larger than the value 0.78 
[22] violate the Clauser-Horne-Shimony-Holt (CHSH) inequality and are thus 
non-local. The entanglement capability C is the smallest eigenvalue of the 
partially transposed density matrix for all possible unentangled input states 
l^m)- ( see below). It has been shown [28] to be negative for an entangled state. 
This quantifier should be -0.5, e.g. for the ideal CNOT, thus characterizing 
a maximally entangled final state. 

In Fig. 5, the deviations due to decoherence of the gate quality factors 
from their ideal values are shown. Similar to most of the rates, all gate quality 
factors saturate at temperatures below a threshold set by the qubit energy 
scales. The deviations grow linearily at higher temperatures until they reach 
their theoretical maximum. Comparing the different coupling scenarios, we 
see that at low temperatures, the purity and fidelity are higher for the case 
of one common bath, but if temperature is increased above this threshold, 
fidelity and purity are approximately equal for both the case of one common 
and two distinct baths. This is related to the fact that in the case of one 
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Fig. 5. Log- log plot of the temperature dependence of the deviations of the four 
gate quantifiers from their ideal values. Here the temperature is varied from w to 
2 • E e . In all cases a — a\ — — 10 . The dotted line indicates the upper bound 
set by the Clauser-Horne-Shimony-Holt inequality 

common bath all relaxation and dephasing rates vanish during the two-qubit- 
step of the CNOT (sec [25] for details), due to the special symmetries of 
the Hamiltonian, when temperature goes to zero as discussed above. Still, 
the quantum degree and the entanglement capability tend towards the same 
value for both the case of one common and two distinct baths. This is due to 
the fact that both quantum degree and entanglement capability are, different 
than fidelity and purity, not defined as mean values but rather characterize 
the "best" possible case of all given input states. 

In the recent work by Thorwart and Hanggi [22], the CNOT gate was 
investigated for a Uy 1 ® cffl coupling scheme and one common bath. They 
find a pronounced degradation of the gate performance with gate quality fac- 
tors only weakly depending on temperature. If we set the dissipation and the 
intrinsic energy scale to the same values as in their work, we also observe 
only a weak decrease of the gate quality factors for both the case of one com- 
mon bath and two distinct baths in the same temperature range discussed 
by Thorwart and Hanggi. However, see Fig. 5, overall we find substantially 
better values. This is due to the fact that for & y ® & y coupling, the Hamilto- 
nian does not commute with the coupling to the bath during the two-qubit 
steps of the pulse sequence, i.e. the symmetries of the coupling to the bath 
and the inter-qubit coupling arc not compatible. The dotted line in Fig. 5 
shows that already at comparedly high temperature, about 20 qubit ener- 
gies, a quantum degree larger than Q w 0.78 can be achieved. Only then, the 
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Clauser-Hornc-Shimony-Holt inequality is violated and non-local correlations 
between the qubits occur as described in [22]. Thus, even under rather mod- 
est requirements on the experimental setup which seem to be feasible with 
present day technology, it appears to be possible to demonstrate nonlocality 
and entanglement between superconducting flux qubits. 

5 Summary 

It has been outlined, how one can model the decoherence of an electromag- 
netic environment inductively coupled to a superconducting flux qubit. We 
have exemplified a procedure based on analyzing the classical friction induced 
by the environment for the specific case of the read-out SQUID. It is shown 
that the SQUID can be effectively decoupled from the qubit if no bias cur- 
rent is applied. The effect of the decoherence on relaxation and dephasing 
rates of single qubits has been discussed as well as the gate performance of 
coupled qubits. We have shown that by carefully engineering the impedance 
and the symmetry of the coupling, one can reach excellent gate quality which 
complies with the demands of quantum computation. 

We would like to thank M. Governale, T. Robinson, and M. Thorwart for 
discussions. FKW and MJS acknowledge support from ARO under contract- 
No. P-43385-PH-QC. 
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